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We determine the quantization condition for the energy levels of two interacting particles in a 
finite box in a "moving frame", i.e. one in which the total momentum of pions is non-zero. This 
condition is valid up to corrections which fall exponentially withe the box size, and holds only 
below the inelastic threshold. It is derived using field theoretic methods, using a generalization 
of previous summation formulae relating sums and integrals over momenta. The result agrees 
with that obtained earlier by Rummakainen and Gottlieb using a relativistic quantum mechanical 
approach. Technically, we expand the finite-volume four-point Green function in terms of the 
infinite-volume Bethe-Salpeter kernel, and determine the position of the poles. The final result 
is written in terms of the two-pion scattering phase shift. Our result can be used to facilitate the 
determination of the scattering phase shift, and can be used to generalize the Lellouch-Luscher 
formula relating finite-volume two-particle matrix elements to those in infinite volume. 
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1. Introduction 

Lattice simulations are necessarily performed in a finite volume. Unlike the exponentially 
small errors in lattice evaluation of hadronic masses or matrix elements with at most a single hadron 
in the external states, the errors when two (or more) hadrons are present in external states the finite- 
volume effects decrease more slowly, as powers of the box size, L, and need to be understood in 
order to obtain physical quantities with good precision. The theory of such effects has been fully 
developed for two particles in their rest frame, i.e. with total momentum P = 0. The spectrum of 
such states was worked out in refs. [||, ^, |]] and the finite volume corrections to the matrix elements 
were obtained in refs. [§, q|. 

In this note, we report on a determination of the spectrum of two-particle states with total 
momentum P ^ 0, which we call a moving frame} Full details are given in ref. [^]. 2 

There are several reasons why this extension is important, e.g. the use of a moving frame 
obviates the need for vacuum subtractions in the isoscalar channel, which is otherwise under poor 
statistical control. Further applications are discussed in refs. [{7|, ||]. 

The energy levels of two-particle system are determined from the position of poles of corre- 
sponding correlation function. By expanding the correlation function in terms of Bethe-Salpeter 
kernel, we can locate the terms which potentially generate power corrections. With the aid of 
summation formulae, we can calculate these corrections up to exponentially small errors. After 
rearranging and resumming we deduce the quantization condition, whose solutions generate the 
location of poles. This condition depends on the scattering phase shift and the box size. Thus, 
we can determine the scattering phase shift from a lattice calculation of spectrum. We sketch the 
derivation of the summation formulae in sec. || and of the quantization condition in sec. |3| 



2. Summation Formulae 

Finite-volume effects for two-hadron states originate from the difference between the sums 
over the discrete momenta in a finite volume and the corresponding integrals over the continuous 
spectrum in infinite volume. To determine the relation between them, we start from the Poisson 
summation formula, 



k 



where the summation on the left-hand-side is over all integer values of n = {n\,n2,rfi), with k = 
{2n/L)n. The approximation is valid when the Fourier transform of g(k,), g(r), is non-singular, 
and is either contained in a finite spatial region or decreases exponentially as \r\ — ► oo. We note 
that functions g(k) with these properties have no singularities for real k, and fall off fast enough at 



| A; | — > oo that the integrals in eq. ( |2.1| ) converge. 



1 A generalisation of the results for the spectrum of such states to a moving frame was proposed some time ago in 
ref. [jj. We address this in sec. ^ 

2 The results have been also obtained in ref. [H] using a different method. 
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However, the finite- volume corrections for two-hadron correlators with non-zero total momen- 
tum are contained in summations of the form 

k 

where we assume that q* 2 is such that there is no term in the sum with k* 2 = \k* | 2 = q* 2 and 
that f(k* ) has the properties discussed above. Note that the summation is over the moving frame 
momenta k = (2%/L)n, with n being a vector of integers but the summand is written in terms of the 
centre-of-mass momenta k* using the Lorentz transformation of eq.(|3.4[). It is also convenient, as 
will be apparent shortly, to pull the Jacobian col / (Ok out of the function f(k* ). 

The manifest singularity at k 2 = q 2 forbids the immediate application of the eq. ( |2.1[ ). In order 
to avoid this difficulty, we expand / in terms of spherical harmonics 3 and subtract from S(q*) a 
function which is specially chosen to cancel the singularity: 

k 1 

^ / i!S^fi!!,* w ,, f) (2 . 3) 

where the Jacobian factor (a£/(Ok corresponds to the change of integration variables from the 
laboratory-frame momenta k to the centre-of-mass frame momenta k*. This relation is valid up 
to terms which are exponentially small in the volume. The exponential factors exp[a(q* 2 — k* 2 )} 
(with a > 0) are included so that the subtraction does not introduce ultraviolet divergences. By 
rearranging terms, the required summation formulae for given spherical harmonics is 4 

SU<f) = S,fi ^j^S P^+f^)clM 2 ) , (2-5) 



where 



i a(q* 2 -k* 2 ) r j3m p a{cf 1 -k* 1 ) 

itf^^E^r^r'^^.r)-*^/^-^. (2.6) 



Although the integrand in eq. ( |2.4[ ) has no pole at k = q, in eq. fl2.5| ) we separate it into two terms 
each of which does have such a pole. For consistency, the two terms need to be regulated in the 
same way and the principal value prescription, denoted is a natural choice. Finally, we gather 
all the spherical harmonics and then the summation formula is 

5(<f) = %-i fim{q * )cLiq * 2) - (2 - 7) 

3 Here, we used \fAnY\ m , which simplifies the normalization for / = 0. 

4 This equality holds up to exponentially small corrections. From now on, this will not be stated explicitly. 
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Figure 1: Diagrammatic expansion of the correlator Cp{E). Propagators are fully dressed and normalized 
to unity on shell. K is the amputated two-particle irreducible four-particle correlation function. The circles 
at the ends represent the operator a, renormalized by two factors of \JZ. 

3. Quantization Condition in Moving Frames 

The two-pion spectrum in finite volume can be extracted from the correlation functions of 
composite operators: 

Cp(t) = (0\ o p (t) ^(0,0)10), Cp{E) = J dte- iEt C p {t). (3.1) 

(with time ordering implicit), where a(x,t) is an interpolating operator for two-pion states and 
op{t) is its spatial Fourier transform. We will determine the quantization condition by identifying 
the position of poles of Cp(E). 

The correlation function Cp can be expressed in terms of the Bethe-Salpeter kernel K through 
the series shown in fig. [I]. Since we choose E to lie below the four-pion threshold, there are no 
intermediate states with four or more pions and the finite-volume effects in K are exponentially 
suppressed [||, 0]. The same is true of the dressed single particle propagators [JTJ] _ The only power- 
law volume corrections arise through the the two pion loops, and we now turn to an analysis of 
these corrections. The generic loop integration/summation appearing in fig. [I] is of the form 5 

l v rdko f(k ,k) 

~ L 3 ^J 2n (k 2 -m 2 + ie)((P-k) 2 -m 2 + ie) K ' 

where k = (ko,k) and P = {E,P) are four vectors and we have left out a factor of i = — 1 from the 
propagators which will be accounted for later. The only properties of / that we need are, first, that 
it has no singularities for real k (which holds given our kinematical constraint on E), and, second, 
that its ultraviolet behaviour is such as to render the integral and sum convergent. 

To simplify the pole structure we first perform the ko integration. We choose to close the 
contour of integration so as to pick up the "particle" contribution from the first pole and the "anti- 
particle" contribution from the second: 

m,k) + f(E + co Pk k) (3 _ 3) 




2co k ( (E - (O k ) 2 - (O 2 ) 2(0 Pk {(E + (0 Pk ) 2 




where co k = V~k 2 + m 2 and (Opk = y (P — k) 2 + m 2 . 

5 Although we phrase our discussion in Minkowski space, we note that this same object may be obtained from the 
Euclidean space correlators calculated in lattice simulations by analytic continuation to imaginary Euclidean energy 
(which is the approach used in ref. [&). 
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For the kinematic region of interest, < E 2 — P 2 < 16m 2 , it is straightforward to show that the 
only singularity in I is the explicit pole in the first term inside the braces in eq. (33), which occurs 
at those values of E for which there is a term in the summation with cot + (Opt = E. Since this 
singularity leads to finite-volume corrections which decrease like powers of the volume, we need 
to examine the first term in more detail. 

With Lorentz transformed variables, 



*fi = Y( k \\ -P«>k), k]_=kj_ and a k * = ^k* 2 + m 2 = y{(O k - J3£||) . (3.4) 

we rewrite the first term (which we call l{) in the form 

, 1 1 r 1 /*(£*) _ .1 l y (Ql /*(**) E* + 2(0* k 
l V> E* V 2fflfc E* - 2co* 'l 3 2E* *s> co k q* 2 -k* 2 4(0? ' 



k 



where f* is the function / rewritten in terms of the centre-of-mass variables. We now apply the 
summation formulae from sec. 0, 



. 1 f cPk* f*(k*) E* + 2co* k i ~ ' 

In order to isolate the finite-volume correction, we replace the principal-value integral in 
eq. (3.6) by the corresponding integral with the Feynman is prescription in the propagator and 
a "delta-function" term: 

J _/_L f**L E * +2( °Z | fftotf) [ Y y r<n*\r P (n^ 

l ~ 2E* J {2nf q* 2 -k* 2 + ie 4co* k + SnE* 2E* 

Note that the "delta-function" term picks out the / = part of f*. Observing that the first term 
in eq. (3/7) is exactly the infinite volume expression for I\ in Minkowski space (after retracing the 
steps in the derivation above), we arrive at the finite-volume correction of loop integrals, 

I=I„+I FV , I F y= j^J^- 4e* t t ftn{<f)cU<t 2 )\ ■ 

Since each loop integral (sum) in fig. [j] contains an infinite volume part and the finite volume 
correction, we can rearrange the expansion according to the number of insertions of the latter. 
Ignoring the zero insertion term, which is irrelevant to the pole structure of the correlation function, 
we can represent the series as in fig. ||, leading to the general result: 

C F /(E) = -A' F A+A' F (iM/2) FA+-= -A 1 F —A. (3.7) 

p y 1 \ I ) l + iMF/2 

Here we have taken into account the factor of i dropped from the loop in the previous section, 
as well as symmetry factors of 1/2 arising from the identical nature of the particles. A and A' 
represent the overlap of a with asymptotic two-pion states with definite angular momentum. F is 
a kinematic factor obtained from Ipy: 



Fhm\h,m z = -^^\hh 8 mim 2 -i-rY, £ —j-cf^q* 2 ) [ <m*Yl m Yt m Y h ^A , (3.8) 

° nr ' \ /=0m=-/ 1 J ' 
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Figure 2: Contributions to the volume dependent part of the a correlator, Cp V . The notation is as in fig. [j] 
except that the filled circles represent the full scattering amplitude, M, given by a geometric sum of any 
number of insertions of the kernel K, and the vertical dashed lines indicate that the on-shell finite volume 
part, Ifv, has been used for the loop integral. The quantities A and A' are defined in the text. 

and M is the on-shell scattering amplitude: 

M XX l67lE * (exp[2ig/,(g*)]-l) 

Mi im; i 2tm2 = 0^i 2 0„ umi — — — , (3.9) 

Now, we can determine the final quantization condition: 

det(l + /MF/2) =0. (3.10) 

Solving this after truncation of the partial wave expansion leads to results which can be shown to 
be equivalent to those of ref. [{7|]. 

4. Conclusion 

We have provided a field theoretic derivation of the finite volume energy shift for two hadron 
states in a moving frame, confirming the result obtained by ref. [|7|]. Our result can be used to 
determine the finite- volume corrections to matrix elements of local composite operators with an 
initial and/or final state consisting of two hadrons, thus generalizes the Lellouch-Liischer factor to 
moving frames [|| |9|]. The path is now open to numerical studies of two pion energies and matrix 
elements in a moving frame, which, as discussed in refs. [0, || §], simplifies a number of 
important applications. 
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